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ABSTRACT 


The purpose of this research was to show, both experiment- 
ally and analytically, that the stress concentration around 
holes in laminated fibrous composites is not adequately de- 
scribed by the theoretical solution for homogeneous orthotropic 
plates, and that it is a function of hole size. It was shown 
that the gross laminate properties, determined analytically from 
individual lamina properties, make it impossible to express tne 
proper boundary conditions at a free edge. A new expression 
for the modulus of elasticity tangential to a free boundary 
was developed. Thus, the problem became similar to that of a 
circular hole strengthened by an elastic ring. Furthermore, 
the effective width of such a ring was found to depend upon 
hole size. Thus, the stress concentration factor was shown to 
be a function of hole size, which was not predicted by homogene- 
ous orthotropic theory. These results were very satisfactorily 


verified by experiment. 
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G = shear modulus 

G/E - glass/epoxy 

E = modulus of elasticity 

K = strain concentration 

n = material parameter given by eqn (II.E.2) 

DON — transverse distance from center of hole (II.E.1) 

Y - longitudinal distance from center of hole (II.E.1) 


€ 


stiffness matrix 


[С] 


[С] 


transformed stiffness matrix 


[T] 


transformation matrix (Eqn III.A. 3) 


€ - strain 
O - stress 
V - Poisson's ratio 
т - shear stress 
Y - shear rem 
а - lamina orientation angle from longitudinal axis 
6 - angle from longitudinal axis 
= load | 
h - thickness 
SUBSCRIPTS 
I = longitudinal axis of composite plate 
D - transverse axis of composite plate 
45 - 45° between 1 and 2 directions 
N - longitudinal (fiber) axis of a lamina 





E - transverse axis of a lamina 


n = normal to plate boundary 

S = tangential to plate boundary 

0 - arbitrary angle from longitudinal axis 
SUPERSCRIPTS 

C - composite plate (includes photoelastic coating) 
с-Р - composite plate less photoelastic coating 

P - photoelastic coating 

(К) - denotes xem lamina 





I. INTRODUCTION 


A. BACKGROUND 

Because of their high strength-to-weight ratios, laminated 
fibrous composites are expected to be utilized extensively for 
aerospace vehicles in the future. At the present time, how- 
ever, the realization of the maximum potential of composites 
is severely restricted by the limited techniques available for 
their analysis. An aerospace vehicle is necessarily a very 
complex structure which is difficult to analyze even when com- 
posed of homogeneous isotropic materials. The introduction of 
inhomogeneous anisotropic composites greatly magnifies the 
problem.  Lekhnitskii presents a solution for stresses around 
elliptic holes in homogeneous orthotropic plates, but the ap- 


plicability to nonhomogeneous composites is questionable. 


Б OBJECTIVES 

The first objective of this research was to show experi- 
mentally that the homogeneous orthotropic solution for stress 
concentration around a circular hole is not accurate for 
laminated fibrous composites. The second objective was to 
determine whether hole size has an effect on the stress 
concentration in an infinite plate. The third objective was 
to analyze the effects of inhomogeneity on stress concentra- 
tion. Specifically, the boundary conditions on a free edge 


in an inhomogeneous laminate were analyzed in general and 


l0 





in particular for the composite tested experimentally in this 


research. 


C. APPROACH 

As stated above, the experimental objectives required . 
the comparison of stress fields for different hole sizes with 
the homogeneous orthotropic solution. A (0,5459), glass- 
epoxy composite was used. Four different hole sizes were 
tested by reflection photoelasticity. The analytical approach 
was to examine the relationships between constituent, lamina, 
and laminate material properties given by Calcote”. The pur- 
pose of this analysis was to determine if this attempt to re- 
present laminated fibrous composites as homogeneous materials 


was valid near free boundaries. 
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il. EXPERIMENTAL ANALYSIS 


A. SPECIMEN DESIGN 
l. Material Properties 

The solution given by Lekhnitskii [3] for the recs 
concentration around a circular hole in an orthotropic plate 
is dependent upon the material properties. Therefore, the 
first task was to determine these properties experimentally. 
Due to the thickness of the photoelastic coating to be used 
for the stress concentration analysis, the entire glass-epoxy- 
photoelastic system was analyzed as a whole, rather than try- 
ing to apply correction factors. Thus, the material properties 
were also obtained by photoelasticity. 

Determination of the material properties of an ortho- 
tropic material requires tension specimens cut along the 
principle 1-Е axes, and a third specimen cut at 45? to these 
axes[2]. Therefore, standard tensile specimens as described 
in the Advanced Composites Design Guide [1] were made up. 

2. Stress Concentration 

Previous experimental work by this author seemed to 
indicate that stress concentration in composites is a func- 
tion of hole size. One purpose of the present research was 
to attempt to verify this hypothesis. Therefore, hole size 
was the only independent variable considered in designing 
specimens. All specimens were cut from the same sheet of 


material. The width of each specimen was twelve times the 
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hole radius. The dimensions of all specimens are given in 


Table II.A.1. 


TABBESTT.2.] 


1200 SERIES SPECIMENS 


MATERIAL:  GLASS/EPOXY/PS1-C 


G/E LAYUP: (0?, л d 


BEHRCKNESS: .15" 


PLATE + LENGTH WIDTH HOLE STRESS ANGLE 

L (in) W (in) DIAMETER (a) 
(2a) 

1200 9 1 = 0° 
1210 9 1 - 45° 
1220 9 ii - 90° 
1203 18 S ү? 0° 
1205 12 1-1/2 1/4 0° 
1206 7-1/2 3/4 1/8 0° 


B. SPECIMEN PREPARATION 

The G/E composite was obtained from 3-M Company in a two 
by three foot sheet. The specimens were cut from this sheet 
using a diamond edge saw. All cutting operations were water 
cooled. The ends of each specimen were built up by attaching 


pieces of (0°, 90°) G/E as described in the Composites Design 
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Guide. The most difficult operation was that of drilling 
holes in the G/E material. The difficulty lay in preventing 
Erber pullout as the drill bit came out the back of the 
material. This problem was overcome by sandwiching the 
composite between aluminum when drilling. 

The PS-1C photoelastic coating was cut one-eighth inch 
oversize and the holes were drilled. After cementing the 
coating to both sides of the G/E composite, the holes were 


reamed and the edges were trimmed by machining at high speed. 


C. TESTING PROCEDURE 

The photoelastic data was taken using the 030 series re- 
flection polariscope by Photolastic, Inc. The specimens 
were mounted in a Reihle 10,000 lb. testing machine. in 
order to reduce the scatter due to viscoelastic effects, 
specimens 1200, 1210, and 1220 were cycled to 300 lb. three 
times before testing began. The other specimens were held 
at their test load for five minutes before taking data. The 


tare values were taken after this initial stretching. 


D. MATERIAL PROPERTIES RESULTS 

The original and reduced data are given in tabular form 
in the Appendix. Note that the fringe value, F, has been 
divided by the scale factor of 47 of the linear compensator. 
Thus, the "fringe orders” given (i.e. NN, NO, TNN, and TNO) 
are actually compensator readings. Using these data, the 
stress-strain curves for each specimen and the Poisson's 


Ratio vs. strain curves for specimens 1200 and 1220 were 
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point plotted (Figures II.D.1-5). There is considerable data 
scatter, particularly in Poisson's Ratio. Although a good deal 
of effort was directed towards eliminating this scatter, visco- 
elastic effects and slight variations in load during testing 
made the task extremely difficult. It appears that the visco- 
elastic properties of the material affect the transverse 

strain to an even greater degree than the longitudinal strain. 
In view of these difficulties, due consideration was given to 


the relationship, 


БМ Б, 


19217 52712 (II.D.1) 


when the curves were drawn in these figures. Using the equa- 


tion 





(EE D2) 


the following material properties were derived from the 


Curves: 
ES = 1.24 MSI = 
1 ў = 
ES = 0.78 MSI = 
2 " Моне“ 233 
6.2 = -51 MSI 


Е. STRESS CONCENTRATION RESULTS 
A point-by-point photoelastic analysis was done for plates 


1203, 1205, and 1206. The points chosen were along the X and 
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FIGURE 11,D,2 - Modulus at 45 Degrees 
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FIGURE 11,D,3 - TRANSVERSE MODULUS 





STRESS 


042 034 


MICROSTRAIN 


о НО UNITS. INGH. 
1 ൮) INCH. 


1) വചി 


18 





EF ICU ET. D,4 - Longitudinal Poisson's Ratio 
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FIGURE ІІ,р,5 - Transverse Poisson's Ratio 
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Y axes as shown in Figure II.E.l. The data are tabulated in 
the Appendix as explained in the previous paragraph. The 
longitudinal strain variation along the X-axis from the hole 
to the edge of the plate is plotted in Figures II.E.2-5. 

The theoretical strain concentration at the hole is given 


by Lekhnitskii as 


К = р(1+п) (ВЕСЕ, 
where 
/ E E 
1 T 
п = 2 (— – у,.) + — (ET Ee 22 


These equations give a theoretical strain concentration factor 
of 3.15 for this material. It is evident from Figure II.D.5 
that the strain concentration is a function of hole size for 
this laminated fibrous composite, even though the hole-size- 
to-plate-width ratio is the same. Also, the measured value 
Seems to approach the theoretical value as hole size increases. 
Previous tests by this author on plates of constant width gave 
the same results. 

Figure II.D.5 also shows that the higher the strain con- 
zeutration at the hole, the faster it drops off with the 
normalized distance from the hole. In fact, the former is 
actually lower at some points! This phenomenon is necessary 
so that the average value of the stress concentration over 
the entire xis (including the hole) will be equal to one. 


111 1196 drop in strain concentration at the outside edge 
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FIGURE 11,E,1 - Coordinate System 








NORMALIZED DISTANCE, x/r 


feet = 2.00700 UNITS INCH. 
tub DUECÜD UNITS INCH. 


PLATE # 1203 «-BXIS 


FIGURE IiI,E.2 ~ Strain Concentration 


SABADO 2 


23 











К 





STRAIN CONCENTRATION, 


NORMALIZED DISTANCE, x/r 


M-SCRLEs2.00€400 UNITS INCH. 
Y-SCALE=S.006-OL UNITS INCH. 


PLATE # 1205 4- 415 


FIGURE ТТ.Е.З - Strain Concentration 


SABA0072 


24 





Е 


К 
25 


23 










STRAIN CONCENTRATION, 


| 


ES 
E 

N 
=> 
съ 
2 


25 | 15 
Т» 
Ca 
г> 7 
съ 
съ 
мъ 


-002 


NORMALIZED DISTANCE, x/r 


В К -2. ВОЕ+00 UNITS INCH. 
ша Eso. 00E=01 UNITS INCH. 


PLATE # 1206 പ 


FIGURE I1,E,4 - STRAIN CONCENTRATION 
BS а m 
SHBHUUT 7 


25 





K 


STRAIN CONCENTRATION, 


ST 


2 


o 


= ИА CSA cc 
ഥ്‌ 


-----|/295 


hb d 1106 


D 2 4 6 


NORMALIZED DISTANCE, x/r 


tits eran Concentration 


26 





of each plate. This effect is due to a change in gross 
laminate properties at a free edge, which is discussed in de- 


tail in Section III. 
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III. THEORETICAL ANALYSIS 


The solution given by Lekhnitskii for the stresses around 
elliptic holes in orthotropic plates assumes a homogeneous 
material. Because the stress distribution is a function of 
the material properties, it is necessary to formulate gross 
laminate properties which approximate the elastic properties 
of the composite by those of a homogeneous system. The obvi- 
ous question here is, “Are the basic assumptions necessary 
for this formulation valid throughout the plate, particularly 


near a free boundary?" 


A. BASIC FORMULATION 

The basic formulation will first be derived in detail be- 
cause it is the basis for the discussion which follows. Тһе 
stress-strain relationship (referred to natural axes) for an 


orthotropic lamina is 


for £C] {Е} (III.A.1) 
where 

{o} = fo ) 

o 1724 

"P 
{e} = Е 

Et 

E 


and 


ot 





44) 


Referred to arbitrary axes, this relationship becomes [2] 


a EE 
с. = [С] => (GILA?) 
A 
12 1/27,, 
where 
= -1 
[C] = [T] (EAS 
and [T] is the matrix of direction cosines: 
2 22 
COS’ a сіп о 2sing cosa 
[Т] = п cO s -2sina cosa (TETAS) 
-sinqa cosa sina Cosa koca M^ 


Equilibrium requires that (assuming unit width) 


и о 


where в (К) is the thickness of the geh lamina. Thus, 
(k), (k) 
T Pn (III.A.4) 
;n 


Substituting equation III.A.2 into this expression yields 


SENE т nD 


c 


(III.A.5) 
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At this point, the assumption is made that the strain field 
in each lamina is identical to that of the composite as a 


whole, that is 


le) = (e)? = CS (СТЕ АБ) 


equals a constant over the summation. Thus, equation III.A.5 
reduces to 


з (к) (К) 


[с] С = E 
zy V) 


(III.A.7) 


B. EREE EDGE BOUNDARY CONDITIONS 
The boundary condition for a free edge in a homogeneous 


plate loaded in-plane is 


o = 0 СЕТИВА 


By applying this boundary condition to a composite represented 

by the homogeneous properties given by III.A.7, one is actually 
only requiring that the average stress across the thickness be 

zero, rather than the stress at each point across the thick- 


ness of the plate. Applying equation III.A.4 to III.B.l, 
COME) 


E Ха h 
o = а = (б (AR BZ) 
д Eh 
Thus, E in general need not be zero, whereas the physics 


of the problem require that it must be zero. Therefore, the 
material properties given by equation III.A.7 are not valid 


near a free boundary in a laminated composite. 
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Equation III.B.2 results in the free-edge displacement 
shown in Figure III.B.2. The free edge is still a plane nor- 
mal to the surface. Thus, all strains and displacements are 


(к) in a lamina 


constant through the thickness. The stress ге 
is derived as follows. The tangential stress, Or is known 


and the normal stress, gu is zero. Thus, at the free edge 


ഞ്ഞ... 
E i ín o / E. (ТТТ.В.3) 
апа 
ко єс | с с S 
€ = € = -у € 
n n com © 
SUE 
a TR ee Wey 
En l-v__V n sn 
ns sn 
or 
E (к) 
(К) _ n ес 
On Я Te Wen oe TTB 
ns sn 


Thus, the error in the boundary condition is directly propor- 
tional to the difference in Poisson's ratios between the 
composite and the lamina in question. The error is zero for 
homogeneous materials. 

The free edge actually looks like Figure III.B.3. The 
displacement must be continuous between lamina. Thus, the 
displacement will vary through the thickness of a given 


lamina. 
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Pigi TTB NE Coordinate System 


ЕТЕТ | 


LB ж Displacement for о” = 


в 


Pig, ва Вон Displacement for stk). 0 


= 


>> 





C. MATERIAL PROPERTIES AT A FREE BOUNDARY 


It has been shown that equation III.B.2 is not the correct 


boundary condition for a free boundary in a laminated compos- 


mee. The correct condition is 


іш — №0 for all К 


n 


(III 


In order to derive the composite material properties, 


tional assumption is required: 


ПЕ У тс 
s 65 T Es 


is constant through the thickness. 


СО Е 
En ба 


Hooke's Law for this case is 


С да 
Os Е. 
and 
(к). a(k) 
A а 


Recalling equation III.A.4 


MER e 
Os = (К) 
OY 
(k) (k) 
n _ LE. eh 
575 yh СК) 
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Note that 


(III 


(IFI; 


(III 


(III 


ae) 


an addi- 


CE) 


с. 


254) 


26 25) 





Thus : 
ae) К) 
S 


tj 
Q 
| 
3 
СЪ 
D^ 


S тъ СК) 


Dee NUMERICAL RESULTS 

The next step in the analysis was to find numerical values 
for ЕС, Young's modulus at a free boundary and i Young's 
modulus in an arbitrary direction far from any boundaries. 
E- may be found by use of equation III.C.6, which requires that 
B" be known for each lamina. A single lamina composite 
could have been made up and tested, but the data scatter in- 
herent in the properties of glass-epoxy composites cut from 
different sheets would have rendered the results virtually 
useless for this analysis. Fortunately, this procedure was 
not necessary. It was possible to analytically reduce the 
gross laminate properties to individual lamina properties by 
the procedure which follows. The experimentally determined 


properties of the G/E/PS1-C composite are repeated here for 


convenience. 


E, = 1.24 MSI 912 = .5 
ой C = 

Е, = ./8 MSI У-1 z o3 
C — 

G32 = .51 MSI 


The stiffness matrix may be computed from 
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E E 


11 1-у12У21 22 1-у12У-1 
ра 
22 
C == — С = С (ТЕТ. БЕ) 
12 1 У12У21 44 12 | 
which give 
1.47 2162 
[C] = .462 .923 0 
0 0 ой ТЕ) 


The properties of the PS1-C photoelastic coating are 


ЕР - .462 
PS 36 
and the stiffness matrix is 
7528 „190 0 
ke =| |150 „528 0 
0 0 „169 (ITE DSS) 
The thickness of the С/Е and PS1-C portions of the composite 


will also be required. 
h` = „15 


um = 0625 


y 
ща 
॥ 


20875 

The thickness of the PS1-C includes the cement, which has ap- 
proximately the same properties. Recalling equation III.A.7, 
it is possible to determine the properties of the G/E system 


without the photoelastic coating. 


Һе [с] = hP[c)P+no7P [cy o7P (III.D.4) 
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yields 


ІСІС?Р = . 843 1.48 0 
0 0 .987 (TTT DEO) 


The inverse of equations III.D.1 are [2] 


2 
E = С - MID 
l ll C22 
2 
C 
12 
E = 16 ren 
2 22 Cii 
С 
Vio = = (ITI 5) 
22 
V = С\2 
21 a 
11 
612 = Opa 


which gives the properties of the G/E system. 


Сари. SM 
Cpu Sap 
С-Р _ 


Recalling equation III.A.7 and noting that the thickness of 


each lamina in the G/E system is the same: 


=. (К) 
[217 = HO a (711010) 


where n is the number of lamina. 
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Recalling that 


= =! 

[С] = [т] Ср Га CETT DIS) 
it is possible to solve for [C] because all lamina are the 
same except for orientation! Equation III.D.8 may also be 


written [2] 


=C Сос (с +2C ЕЕ SS 


പി io any 22 

С = (CTC AAC es oo ole er cco so) 
12 Jl 22 4 4 12 

C=C C +2C ee sens oa (IT-D. 9S) 
РЕ : 12 44 22 ‘ а; 

€ =(C,,+C__-2C, „-2C ന്‌ ദോ IG O 
44 O 12 44 44 


Noting that the C15 5 do not vary across the summation in 
equation III.D.7, one may proceed as follows. The orientation 


Of the various lamina are 


ШЕ) (6) 


а = а = 0° 
gU? = a = +45° 
ШЕ) ЖШС) е 


and equation III.D.7 reduces to 





е-р Е | 
Cii 3 2 1 Cii 
Ci» l 4 l С,» 
= с (III.D.8) 
C22 1 2 3 C22 
wr WEE з Сад 
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The inverse is 


С iS 1 
om „5 0 
Cae laos 1 
Сад | ട്‌ 


which gives the stiffness 


4.78 
163 
0 0 


s161 


[C] = 7952 


с-р 
2 Cii 
-2 Су. 
(тт. р.9) 
2 C22 
al Los 


matrix for a single lamina as 


(III.D.10) 
‚ 305 


Applying equation III.D.6 gives the material properties of 


the lamina. 


E€ = 4.75 MSI 
E - .847 MSI 
Got” 2505 MSI 


In order to apply equation III.C.6, E. 


each lamina. 





gu € 
e Е eee: + ar sums 

Е (К) a 

S 
Or 
55 = .211cos^?? + 1.18sin 
УЕ 
S 


"E 


tz 
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d 
-- 1/4 [> = 27 1Sin 2 
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. 03 


(К) must be known for 


This quantity may be computed from 
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FIGURE III,D.l - Young's Modulus Curves 
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Е is computed for each lamina and then ЕС is computed for 


equation III.C.6. 


For comparison, ES may be computed from 








1.44 = TEN. Ө + y а 9 + „358 sm (20) 
E 
0 
or 
1 = .806 соѕ Ө + 1.28 зіп? Ө + .289 sin” (20) 
Е (ED ON 


5 and E- are plotted as a functies ©f 0 in Figure III.D.l. 


E. EFFECT ON STRESS CONCENTRATION 

As seen in Figure III.D.1, Young's modulus in the y direc- 
tion at the boundary is 1.17 MSI compared to 1.24 MSI far from 
all boundaries. The variation along the X-axis is something 
like Figure III.E.1. Thus, the material at the boundaries is 
less stiff than the material near the center of the plate. 
Now, suppose that the theoretical stress concentration varia- 
tion for constant material properties is given by curve A in 
Figure III.E.2. The variation in material properties described 
by Figure III.E.1 would tend to reduce the stress concentra- 
tion ERE the boundaries, because the material is less stiff 
in these areas (Curve B). This effect is shown clearly at 
the boundaries of the plates in Figures II.D.2-5. Thus, the 


lower stiffness near a free boundary in a laminated composite 
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requires that the stress concentration factor be smaller than 
that given by homogeneous orthotropic theory. 

The effect described above can also be shown to vary with 
the hole size. The modulus of elasticity of the laminated 
composite is lowest at the free boundary and increases to a 
steady value at a finite distance from the boundary as shown 
in Figure III.E.l. The actual function will depend on compon- 
ent materials and layup of composite, but will depend only on 


the X distance from the free edge for a particular composite: 


Е = (X) TEE SET; 


The stress concentration along the X-axis, however, is a 


function of the number of hole radii away from the boundary: 


K = K(X/r) (TET 2) 


Thus, in terms of the normalized distance, X/r, a plate with 
a small hole is affected more by the difference in stiffness 
mona plate with a large hole. Therefore, the smaller the 
hole, the smaller the stress concentration factor, even in an 
infinite plate or plates with a constant hole-size-to-plate- 
width ratio. 

If Young's modulus in Figure III.E.l is assumed to vary 
like Curve B, rather than Curve A, one has the problem of 
stress concentration around holes which are strengthened by 
elastic rings. This problem has been solved by Savin [4] who 
describes similar results. The theoretical stress concentra- 


tion for infinite copper plates with holes strengthened by 
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R free boundary X 


FIGURE 111,E,1 - E near a free boundary 


——— 
X 


FIGURE 111,E,2 - Inhomogeneous effects on k 
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steel rings of various widths is presented by Savin and this 
set of curves is very Similar to Figure II.E.5, the experi- 
mental results for plates with different size holes. Consider 
a series of such rings, one inside the other. A large number 
of thin rings with appropriate material properties could be 
used to produce the variation in stiffness described by Curve 
A of Figure III.E.1. It is suggested that this approach be 


further investigated in order to obtain numerical results. 
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IV. CONCLUSION 


The first objective of this research was to show experi- 
mentally that the homogeneous orthotropic solution for stress 
concentration around a circular hole is not accurate for 
laminated fibrous composites. Referring to Figure III.E.5, 


the following stress concentration factors are obtained: 


TABLE IV.l 


Theory о = 
Plate 1203 (r = 25") 3.07 


Plate 1205 (r 


125.) 2.45 


Plate 1206 (r 


0 101 


Obviously, the theory is not accurate. The second objective 
was to determine whether hole size has an effect.on the stress 
concentration in an infinite plate. Referring to Table IV.1 
again, it is obvious that the stress concentration factor de- 
creases with hole size. 

The third objective was to analyze the effects of in- 
homogeneity on stress concentration. It was found that to en- 
force the boundary condition at a free edge required the 
Меса топ of the gross laminate properties near the edge. 
The lower stiffness obtained analytically required that the 
stress concentration factor be lower than that predicted by 
homogeneous orthotropic theory. Furthermore, it was shown 


that the stress concentration factor varies with hole size 
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because the modified material properties affect a larger 


normalized area for smaller holes. 
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